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^ ■ ABSTRACT 

o 

I We report detailed studies of the X-ray emission from the kiloparsec-scale jet in 



the nearest active galaxy, Centaurus A. By analyzing the highest quality X-ray data 



Qh' obtained with the Chandra ACIS-S, 41 compact sources (mostly bright jet-knots) were 

O ' found within the jet on angular scales less than 4 arcsec, 13 of which were newly identi- 

' fied. We construct the luminosity function for the detected jet-knots and argue that the 

■ remaining emission is most likely to be truly diffuse, rather than resulting from the sum 

of many unresolved fainter knots. We subtracted the contributions of the bright knots 
^ , from the total X-ray jet flux, and show that the remaining extended emission has a 

H ' relatively flat-topped intensity profile in the transverse jet direction, with the intensity 

peaking at the jet boundaries between 50" and 170". We note that limb-brightened 
morphologies have been observed previously at radio frequencies in a few FR I and 
FR II jet sources, but never so clearly at higher photon energies. Our result therefore 
supports a stratified jet model, consisting of a relativistic outflow including a bound- 
ary layer with a velocity shear. In addition, we found that the X-ray spectrum of the 
diffuse component is almost uniform across and along the jet, with an X-ray energy 
spectral index of qx ~ 1, similar to those observed in the compact knots. We discuss 
this spectral behavior within a framework of shock and stochastic particle acceleration 
processes, connected with the turbulent, supersonic, and non-steady nature of the rel- 
ativistic outflow. We note some evidence for a possible spectral hardening at the outer 
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sheath of the jet, and manifesting itself in observed X-ray spectra of ax < 0.5 in the 

most extreme cases. Due to the hmited photon statistics of the present data, further 
deep observations of Centaurus A are required to determine the reality of this finding, 
however we note that the existence of the hard X-ray features at outer jet boundaries 
would provide an important challenge to theories for the evolution of ultra-relativistic 
particles within extragalactic jets. 

Subject headings: acceleration of particles — galaxies: jets — galaxies: individual 
(Centaurus A) — X-rays: observations — radiation mechanisms: nonthermal 

1. Introduction 

The Chandra X-ray Observatory has confirmed that X-ray emission from large-scale jets is 
common in radio galaxies and quasars (e.g., Harris &: Krawczynski 2002; Kataoka & Stawarz 
2005, and references therein). In nearby FR I sources, the typical radio-to-X-ray spectrum of 
kpc-scale jet knots is consistent with a single, smoothly broken power-law continuum, suggesting 
that this broadband emission is entirely due to nonthermal synchrotron radiation from a single 
electron population (e.g., Hardcastle et al. 2001 for 3C66B), although spatial ofi^scts observed in 
these sources between the positions of the X-ray and radio knots may indicate a more complicated 
physical picture. Meanwhile, the observed X-ray knots in hundred-kpc-scale quasar jets are much 
brighter than expected from a simple extrapolation of the radio-to-optical fluxes (e.g., Schwartz et 
al. 2000 for PKS 0637-752). This has led to the hypothesis that the X-ray production involved 
inverse-Comptonization of CMB photons within a highly rclativistic flow (Tavecchio et al. 2000, 
Celotti et al. 2001). It has been argued, however, based on recent studies of multi- wavelength 
properties of powerful large-scale jets, that the inverse-Compton model may face serious problems 
(see, e.g., the discussion in Stawarz et al. 2004; Kataoka &: Stawarz 2005). Moreover, the inverse- 
Compton model cannot readily explain the detected X-ray emission from a few FR II jets — 
Pictor A (Wilson et al. 2001), 3C303 (Kataoka et al. 2003), 3C 120 (Harris et al. 2004), 3C403 
(Kraft et al. 2005) — since they are expected to be the de-beamed analogues of the radio-loud 
quasars. Therefore, the synchrotron origin of the keV photons remains a plausible alternative even 
for powerful FR Il/quasar jets. 

The synchrotron model for X-ray emission from the large-scale jets in radio galaxies (FR I/FR II) 
and quasars involves much lower bulk velocities for the fiows, but on the other hand requires the 
presence of ultra-relativistic electrons with energies Eg > 100 TeV. This requirement is supported 
in a natural way by the expectation that the large scale jets are good candidates for the the accel- 
eration of cosmic ray particles owing to their large, extended structures and turbulent nature (see 
Hillas 1984). From the theoretical perspective, it has been proposed that the most energetic parti- 
cles are accelerated preferentially in the turbulent boundary layers of such outfiows (Ostrowski 2000; 
Stawarz &; Ostrowski 2002; see also in this context Rieger &; Duflcy 2004 and references therein). 
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In fact, the 'limb-brightened' jet morphology expected in this type of models has been observed at 
radio frequencies (probing relatively low electron energies) in several objects, such as M 87 (FR I, 
Owen et al. 1989) and 3C353 (FR II, Swain et al. 1998). More importantly, the continuous acceler- 
ation processes acting thereby should result in the formation of flat spectra at the highest electron 
energies, which may account for the hard X-ray spectra and high X-ray luminosities of powerful 
jet sources. Similarly, Dermer & Atoyan (2002) suggested the formation of high-energy spectral 
hardening owing to inefficient Klein-Nishina radiative losses of electrons, which also could appear 
at X-ray energies for some particular choices of jet parameters. Therefore the investigations of (1) 
the transverse structure of X-ray jets, and (2) the spectral properties of diffuse X-ray jet emission, 
provide important clues for an understanding of the production of cosmic ray particles and the 
origin of hard X-ray emission in large scale jets. Unfortunately, only few such studies have been 
reported so far, mainly due to the technical difficulties. In fact, it is generally impossible to resolve 
jets in low power FR I radio galaxies or distant quasars, even with the excellent spatial resolution 
(0.5" half power diameter on-axis) and sensitivity of Chandra (see Proposers' Observatory Guide 
v.7^). 

Centaurus A (hereafter Ccn A) is the nearest active galaxy (d = 3.4 Mpc, for which 1" cor- 
responds to 16 pc) and has been well studied across the entire electromagnetic spectrum at the 
highest linear spatial resolution (Israel 1998). It is considered to be a prototypical FR I radio 
galaxy, although the one-sidedness of its kpc-scale jet and the presence of the outer, Mpc-scale 
radio lobes is unusual for this type of object. Previous observations of Cen A have resulted in the 
detection of a bright X-ray and radio jet extending > 4' northeast of the nucleus (Schreier et al. 
1979, 1981; Clarke et al. 1986, 1992; Kraft et al. 2002; Hardcastle et al. 2003). From the analysis 
of the broadband radio-to-X-ray jet emission, the synchrotron origin of keV photons has been 
established for the bright knots, which are characterized by a minimum pressure ~ 10-9 dyn cm , 
equipartition magnetic field ~ 100 /xG and 1 — 10 keV luminosities of the order of ~ 10^^ erg s~^. In 
addition, a very complex morphology of the whole outflow was noted in both the radio and X-ray 
bands, including the presence of bright filaments and diffuse extended regions with complex sub- 
structure. The diffuse X-ray emission clearly extends to the edge of the radio jet and is reasonably 
well matched to the radio emission on these scales. Moreover, Hardcastle et al. (2003) noted two 
conspicuous regions where the extended X-ray emission appeared to be associated with the edges 
of the jet, and to be absent in the center. These X-ray edges of the jet have no obvious radio 
counterparts, although Clarke et al. (1986, 1992) also reported a limb-brightened radio morphology 
in some sections of the jet. We note that these radio polarization studies indicated a projected 
magnetic field parallel to the jet axis along the whole jet length. 

Stimulated by these previous reports, the goal of this work is twofold, namely determining 
the transverse profile of the diffuse X-ray emission extending over the Cen A jet, and constraining 
its spectral properties. Our approach differs from previous studies in that we focus on the diffuse 



http: //cxc .haxvard. edu/proposer/POG/ index .html 



-4- 



emission rather than the bright jet knots, which are thought to be compact shocks inside the jet 
(Hardcastle et al. 2003). We remove the contributions from the knots, as particle acceleration 
may take place in quite different manners in different jet regions (but we will revisit this problem 
later in the discussion section). The paper is organized as follows: In §2, we describe the Chandra 
observations, data reduction process and our analysis method. In §3, we present the results of 
the data analysis. In §4, we discuss our findings in the context of particle acceleration models, 
presenting our conclusions in §5. 

2. Data Reduction 

Since the launch of Chandra X-ray Observatory in July 1999, X-ray observations of Cen A have 
been conducted 10 times, three of which were for calibration (CAL) purposes, and the remaining for 
General Observer (GO) or Guaranteed Time Observing (GTO) proposals. The first Chandra X-ray 
image of Cen A was obtained with the microchannel plate X-ray imaging detector HRC-I (Obs 
ID 463, 1253, and 1412), and the results presented in Kraft et al. (2000). The first useful spectral 
information was subsequently presented based on observations using the ACIS-I configuration (Obs 
ID 316 and 962: Kraft et al. 2002). The exposure times for the two observations were 35.9 ksec and 
36.5 ksec, respectively. Although the jet was 4—5' off-axis from the aim point in these observations, 
and the image resolution was therefore significantly degraded, 31 knots in the Cen A jet were 
successfully detected. 

More recent Chandra observations of Cen A were made in 2002 September (Obs ID 2978) 
and 2003 September (Obs ID 3965) using the standard ACIS-S configuration, where the back- 
illuminated S3 chip was focused to obtain maximal sensitivity to soft photons. The exposure times 
for the two observations were 45.2 ksec and 50.2 ksec, respectively. The active nucleus was at the 
aim point, resulting in both sub-arcsecond resolution and good spectral sensitivity for the whole 
inner part of the jet in Cen A. Although the imaging degrades substantially for sources more than 
1' off-axis, and the X-ray jet in Cen A extends over ~4' in length, the encircled energy for a point 
source is contained within a diameter of ~3" all along the outflow. The roll angle of the satellite 
was carefully chosen to avoid the readout streak of the bright nucleus impinging on the jet image. 
The results from the first observation, Obs ID 2978, are given in Hardcastle et al. (2003). The X-ray 
image was compared with new, high dynamic range VLA images for one-to-one identification of 
radio knots, for which apparent sub-luniinal motions were found (on scales of hundreds of pc). Here 
we analyze the combined ACIS-S data (Obs ID 2978 and 3965) to obtain the best image resolution 
and spectral information reported to date. 

A co-added, exposure-corrected image of the Cen A jet in the 0.4— 8keV bandpass is shown in 
Figure 1 (upper). The image is smoothed with a (T = 0.5arcsec Gaussian function. The nucleus is 
located in the southwest corner of the image, and the jet extends to the northeast. The figure shows 
that the diffuse emission from the jet is observed up to ~ 4' from the nucleus. Also, the Chandra 
image shows that the X-ray jet consists of many knots of enhanced emission embedded within the 
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diffuse regions, as described in detail by Kraft et al. (2002) and Hardcastle et al. (2003). Since the 
primary purpose of this paper is to study the characteristics of the diffuse X-ray emission associated 
with the jet, our approach is complementary to previous reports in the literature: we first resolve 
the point-like sources within the jet (in projection). ^ Next we subtract the contributions from 
these point sources to enable a quantitative measurement of the diffuse X-ray emission. One may 
ask if the remaining emission is truly diffuse in nature, or whether it is the summation of numerous 
unresolved jet-knots at even smaller spatial scales. We will revisit this problem in §3.3, based on 
the luminosity function of the jet-knots. 

Using a wavelet decomposition source-detection algorithm {wavdetect; Vikhlinin et al. 1995) 
to detect emission enhancements on angular scales of 0.5" and 1", we find a total of 35 distinct 
knots or enhancements of X-ray emission embedded within the diffuse regions. We then extend 
the angular scales to 2" and 4", but only an additional six sources are found in these larger region 
sizes. This suggests that most of the jet knots in Cen A are very compact and that their extent 
is smaller than, or comparable to, that expected from the Chandra point spread function. Our 
result is therefore consistent with findings by Kraft et al. (2000) and Hardcastle et al. (2003), who 
claimed that most of the X-ray knots are diffuse but at the 1" (~10pc) level. Due to the good 
photon statistics and excellent angular resolution during the observations, the total number of the 
detected sources is significantly larger than that reported in the literature (i.e., 41 compared with 
31 for Kraft et al. 2002). Three of the weaker knots reported in Kraft et al. (2002), NXl (Knot-1), 
EX3 (Knot-19), and FX4 (Knot-23), were not detected at the previous level in our high resolution, 
co-added ACIS-S image. The location, distance and position angle of the 13 knots or enhancements 
which were newly detected in our analysis are listed in Table 1. 

We next apply dmfilth implemented in the CIAO 3.2 package to replace the pixel values for the 

point sources detected by wavdetect with values interpolated from surrounding background regions 
(DIST option in dmfilth; see Science Threads for CIAO 3.2"'). Here the source regions are defined 
as the ellipses output from wavdetect which contain at least a So" excess of photon counts over 
the background. The background is taken from the "doughnut" shaped region surrounding the 
source, with the dimensions of the major and minor axes of the outer ellipse both doubled. The 
resultant image is almost flat, as expected, but contains a few irregTilar pixels with outstanding 
photon counts. This is because, in some regions, source ellipses overlap. The background region for 
a certain point source therefore inevitably contains photons from neighboring source regions. We 
therefore repeat the same process (i.e., source detection by wavdetect and smoothing by dmfilth) 
three times, until no irregular bins appeared in the image. The resultant X-ray image of the diffuse 



^Prom a comparison of the source density within the jet and that in surrounding regions, it is highly Ukely that 

most of sources are jet knots, with only a few likely to be background, or foreground, sources (mainly X-ray binaries 
in Cen A). A check of the NASA/IPAC Extragalactic Database (NED) yielded only one possible identification, with 
an X-ray star. Similarly the chance probability that one of the new X-ray point sources discussed here is a background 
AGN is less than 1% (Kraft et al. 2001). See Table 1 for more details. 

^http : / / cxc . harvEird . edu/ ciao/ threads/ index . html 
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emission is given in Figure 1 {bottom: 0.4— S.OkeV). We also made bandpass images for soft X-rays 
(0.4— 1.5keV) and hard X-rays (1.5— S.OkeV) using the same procedure. 

3. Data Analysis and Results 

In order to derive the intensity (photon counts) profile of the jet both in the longitudinal and 
transverse directions, we first define the position angle (0) and the distance (r) with respect to the 
nucleus of Cen A as shown in Figure 1 (bottom). Note that the position angle is measured from 
east through north (i.e., in a clockwise direction in the Figure). Then we can integrate the photon 
counts in a sector of an arbitrary range of r and 0, to construct the intensity profile parallel or 
perpendicular to the jet main axis, which is given by Q ~ 35°. Since our approach is simply to sum 
the counts in each sector bounded by the stated values of r and 9, and not to derive the emissivity, 
no normalization by the area of each sector (oc 1/r^) was performed. 

3.1. Longitudinal Jet Profiles 

We first integrate the photon counts in each thin shell within a radial distance r ~ r+Ar over 
the range of position angle 25° < < 45° , and then construct the longitudinal intensity profile in 
Ar = 4" steps, from r = 0" to 250". Figure 2 shows the longitudinal jet profiles thus produced, for 
both soft X-rays (0.4— 1.5keV: top) and hard X-rays (1.5— S.OkeV: bottom). The dotted histograms 
show the profiles derived from the raw X-ray images, whereas those for the diffuse X-ray emission 
are given as the data points on the solid-line histogram. We also show the intensity profile of 
the background (i.e., off-jet) region as a dashed line for comparison, where the background was 
accumulated in the 15° < < 25° and 45° < < 55° regions. The intensity profiles along the jet 
are quite different in the two energy bands, especially at the innermost part of the jet (0< r < 50"), 
where the soft X-ray photons are significantly absorbed in the dust lane crossing the Cen A host 
galaxy. 

To see this more clearly, we divided the 1.5— S.OkeV longitudinal jet profile by that of the 
0.4— l.SkcV, to check for hardness variations along the X-ray jet. Figure 3 compares the hardness 
ratio (HR) and the 0.4— S.OkeV X-ray intensity profile. One can clearly see that HR changes 
significantly at rc — 50", corresponding to a linear scale of ~ 800 pc. Obviously, such a dramatic 
change of the HR is not due to the changes of intrinsic jet spectrum. In fact, we can see that the 
HR changes in a similar way for the background region also, as shown by the dashed line in the 
bottom panel. The inner 1 kpc of the jet is observed through the dust lane of the galaxy, making 
detection of any optical component difficult (Israel 1998). This same dust lane, or equatorial gas, 
heavily absorbs soft X-rays as well, resulting in the significant changes of the hardness ratio seen 
in Figure 3. In the outer part of the jet (50" < r < 240"), the HR is almost constant. This trend 
suggests that the X-ray spectrum does not change significantly along the jet, but the contribution 
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from the background photons must be carefully examined before reaching this conclusion. 

To find possible spectral changes along the jet, we fit the diffuse emission with a power- 
law function in XSPEC, for each shell in Ar = 10" steps between 10" and 100", and 20" steps 
between 100" and 240". We do not examine the innermost part of the jet (r < 10") because the 
X-ray emission from the nucleus significantly contaminates the jet emission there. ^ We take the 
background from neighboring shells just outside the jet (15° <Q < 25°, 45° <@< 55° at the same 
range of r) and subtract it from the source spectrum. The result of spectral fitting along the jet 
is shown in Figure 4. As expected from the hardness profile, the X-ray spectrum clearly changes 
at a characteristic distance of Tc 50". In the inner jet region (r < Vc), the spectrum is well 
represented by a power-law function heavily absorbed with a free zero-redshift absorbing column. 
A combined spectrum between r = 10" and 50" gives the best-fit absorbing column density of iVn = 
(4.791q'53)x10^^ cm~^ and an intrinsic energy spectral index of ax = O.SO^'lQ'Jg, where the reduced 
is 0.50 for 164 degrees of freedom (dof). 

The absorption column in the outer region (r > rc) is much smaller, and similar to the Galactic 
column density (~0.7xl0^^ cm~^). Although there is a weak trend of spectral steepening at larger 
distances from the nucleus, there is no sudden change in the spectral form over the wide distance 
of r = 50" and 240". The subtle anti-correlation between the fitted X-ray spectral index and the 
X-ray flux which can be seen in Figure 4 is most likely due to minor variation of the fitted absorbing 
column density at r > Tc. The best-fit parameters for the average spectrum in this outer jet region 
are A^h = (0.96±0.13)xl02^ cm'^ and ax = 1.19±0.07, which have a reduced of 0.88 for 267 
dof. This energy spectral index of the diffuse emission inside the jet is consistent with that reported 
in Hardcastle et al. (2003), viz, A^h = (l.l±0.3)xl02i cm-^ and ax = l-OOl:^:^^, derived for the 
extended jet emission in the jet between knot A2 and the knot B region (approximately 20" < r < 
70"), excluding all compact X-ray features. Interestingly, in the lOOkpc-long jet in the radio galaxy 
Pictor A, the observed X-ray spectral index is constant along the jet, and is ax ~ 1 (Wilson et al. 
2001), as in the kpc-scale jet analyzed here, suggesting the action of similar particle acceleration 
processes for the high-energy electrons resulting in the synchrotron X-ray emission. 

3.2. Transverse Jet Profiles 

We next integrate the photon counts in each narrow sector with an angular width A@ over 
certain radial distances (selected in the range 50"<r<250"), and then construct the transverse 
intensity profiles in AQ = 1° steps, from B = 20° to 50°. We do not examine the innermost part of 
the jet (r < 50") because of the strong absorption due to the dust lane, as discussed in the previous 



''The contamination level is quite uncertain, as the bright nucleus emission is strongly piled up. We however make 
a crude estimate that more than 80% of the photons come from the nucleus, rather than from the jet, within the 
most inner part (r < 10") of the jet. See Fig. 2 for a comparison of background profile with the diffuse jet emission 
profile. 
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section. Moreover, we cannot determine the position angle @ in the inner-jet region with sufficient 
accuracy. For example, the ~ 1" diameter of a typical Chandra image resolution corresponds to 
AG > 1° at r<50", which is larger than the angular step size of AO = 1° for our analysis. Note 
that AO = 1° at r = 250" corresponds to ~4.4", or to a linear scale of ~ 70 pc. Figure 5 shows 
the transverse jet profiles thus produced for both soft X-rays (0.4— 1.5keV: top) and hard X-rays 
(1.5—8.0 keV: bottom). We made separate plots of the transverse profiles for r < 170" and r > 170", 
because the jet's emission decreases to the backgroTuid level at 170", but then gradually increases 
again to a peak at r ~ 210" (see Figures 1 and 2). The dotted histograms show the profiles derived 
from the raw X-ray images, whereas those for the diffuse X-ray emission are given as the data 
points on the solid- line histogram. The transverse jet profile of the diffuse emission shows a flat-top 
structure for the jet region between 50" and 170" but, in contrast, a single, narrow peak for the 
outermost part of the jet (170" < r < 250"). Such an edge-brightened structure, as observed for 
the inner jet region, is not readily apparent in the asymmetric intensity profile derived from the 
raw data. 

The flat-topped intensity proflle along the transverse jet direction in Cen A cannot be explained 
by a homogeneous plasma fluid moving with a uniform velocity. One natural and simple possibility 
to account for the observed profile is to consider a stratified model for the kpc-scale Cen A jet, with 
a velocity shear across the jet outflow. In general, such a structure is actually required to explain 
many observed properties of extragalactic jets, and is also expected on the grounds of theoretical 
and numerical jet studies (see the discussion in Stawarz & Ostrowski 2002 and references therein) . 
In this model, different Dopplcr enhancements of the emission produced within different velocity 
components can result in the observed limb-brightening of the jet. The other possibility discussed 
in the literature is to consider a mildly-relativistic jet bulk velocity and the enhanced emissivity 
of the non-thermal plasma at the jet boundaries resulting from, e.g., non-uniform radial structure 
of the jet magnetic field due to jet interaction with surrounding medium or some particular large- 
scale magnetic topology (see, e.g., the discussion in Owen et al. 1989). Because of the apparent 
one-sidcdncss of the kpc-scale jet in Cen A, we favor the former possibility as the dominant factor 
shaping the transverse intensity profile, although both Doppler favoritism and plasma emissivity 
stratification may both play a role. 

The exact velocity profile across the kpc-scale jet in Cen A cannot be precisely modeled with 
the present data. Let us however underline a few appropriate constraints. First, note that the 
observed parallel magnetic field structure along the jet axis, with no evidence for the transverse 
component that is usually observed in FR I sources, is consistent with a well-developed velocity 
shear (see in this context Laing &; Bridle 2002). Second, if the limb-brightening of the diffuse X-ray 
emission is predominantly due to varying Doppler enhancements of different jet axisymmctric layers, 
and the emissivity is uniform along the jet (although this may not be a realistic assumption for a 
stratified jet), a component for which the emission is most strongly beamed toward the observer is 
the one characterized by the bulk velocity /3(0) = cos a (or bulk Lorentz factor r(0) = 1/ sina), 
where a is the jet viewing angle. With the preferred a ~ 50° (Tingay et al. 1998; Kraft et al. 2003), 
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one can roughly estimate the bulk velocity of the limb-brightened regions (0 ~ 30° and ~ 40°) 
as /3 ~ 0.65 (or T ~ 1.3). This very rough, illustrative estimate indicates that the spine of the 
kpc-scalc jet in Ccn A — at the position angle O = 35° — is most likely still relativistic, with a 
bulk Lorentz factor T > 2. In this context, the single narrow peak at the outermost part of the jet 
(170" < r < 250"), may correspond to the spine, decelerated significantly at the terminal part of 
the coUimated outflow. Finally we note that the half-opening angle of the X-ray jet, 0x < 5°, is 
smaller than that of the radio jet structure (/>r ~ 10° for r > 60" (Clarke ct al. 1992; see also Kraft 
ct al. 2002). In the framework of the above discussion, this would be consistent with the boundary 
shear layer spreading in the transverse direction as the jet propagates. 

Similar to the analysis in the longitudinal direction, we divided the 1.5— S.OkeV transverse jet 
profile by that of the 0.4— l.SkcV, to examine the hardness ratio (HR) variation across the X-ray 
jet. Figure 6 compares the HR and the 0.4— S.OkeV X-ray intensity profiles. Here wc divided the 
inner jet region further into three separate sectors, (a) 50" < r < 80", (b) 80" < r < 120", and 
(c) 120" < r < 170", for a detailed comparison of changes in transverse profile. Edge-brightened 
structTirc is most clearly seen in the profile derived for 50" < r < 80" (Figure 6a) . Interestingly, one 
can clearly sec that the HR is almost uniform across the jet, but shows a weak sign of enhancement 
at the very edges of the X-ray jet regions (Q 25° and Q 45°), followed by a sudden "drop" 
(softening) close to the edges. Similar hardening or softening may explain some variation of the 
HR even for the jet-edges in Figure 6 (b) and (c). We may speculate that the spectrum changes 
considerably at the outer edges of the jet, accounting for these variations, though the present data 
do not allow this to be established. 

First, since the HR contains a significant contribution from background photons, apparent 
changes in HR might be due to variation in the background characteristics rather than changes in 
the source spectrum itself. This effect is important especially at the jet edges, where the emission 
is dominated by background photons. We will carefully examine this problem later by choosing 
an appropriate background for spectral fitting. Second, it is well known that the response of the 
Chandra X-ray mirror changes with the incident X-ray energy. The half-power diameter of the 
point spread function is 0.5" at IkeV, but slightly worse at higher energies (e.g., 1" at 5keV; see 
The Proposers' Observatory Guide. Such an energy dependence of the angular resolution might 
introduce HR variations perpendicular to the Cen A jet, especially at the jet edges. To check 
for this effect quantitatively, we smoothed the 0.4—8.0 keV image (Figure 1 bottom) with Gaussian 
functions of (j = 0.5" and 1.0" to mimic the images obtained in different energy bands. The resultant 
intensity profiles have only negligible (less than 2%) changes in the HR across the jet, suggesting 
that we cannot explain the observed hard X-ray features shown in Figure 6 by this effect alone. 

Finally, variations of HR may be caused by any possible artifact caused by dmfilth and inclusion 
of background. We therefore remove all the point sources detected by wavdetect and fit the diffuse 

emission with a power-law function in XSPEC, without filling the 'holes' remaining after subtracting 
the source ellipses. We assume a fixed absorption column density of = 0.96x10^^ cm~^, which 
is valid all along the jet if r > 50" (see §3.1). We take the background from two sectors just 
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outside the jet (19°<e<23°, 47° < 9 < 51° within an equal r ~ r+Ar) and subtracted them 
from the source spectrum. The result of spectral fitting across the jet with AQ = 2° steps is shown 
in Figure 7. Although the photon counts are inevitably poor at the jet edges, the general trends 
of the flux variations and spectral changes are very similar to that expected from the HR profiles. 
The spectral index seems to increase slightly and gradually along the main jet, being ax > 1-0. 
In addition, the hardest X-ray spectra are observed at the edges of the jet in the 50" < r < 
80" region, where ax = 0.231^;^? for e~25° or ax = 0.40tg|^ for e~45°. We note that the 
background spectrum is approximately represented by a power-law with an energy spectral index 
of ax ~ 1.5—2.0, and therefore background contamination cannot produce an X-ray spectrum as 
hard as that as observed in Figure 7(a). Unfortunately, the photon statistics are insufficient to 
confirm that the X-ray spectrum is actually harder at the very edges of the jet. Nevertheless, this 
provides an important motivation for future deep observations of the Cen A jet. 

3.3. Diffuse or Unresolved? 

As pointed out by Kraft et al. (2002), the X-ray morphology of the Cen A jet is much more 
complex when observed at higher spatial resolution. This raises the question whether some fraction 
of the extended emission discussed in this paper can be explained by the pile-up of small scale 
knots below the resolution limit for Chandra point source detection. To estimate what part of the 
remaining unresolved emission of the jet could indeed result from a large number of unresolved 
low luminosity compact features, we construct a luminosity function (LF) for the already resolved 
jet-knots, presented in Figure 8. Here we consider the 32 jet-knots located between 50" and 250" 
as the strong absorption due to the dust lane makes determination of the unabsorbed luminosity 
quite uncertain for r < 50". We assume a power-law energy spectrum with the spectral index ax 
= 1.0 (Hardcastle et al. 2003; Kraft et al. 2002) and A^h = O.OGxlO^^ cm^^ (see §3.1) to estimate 
the 0.5-5 keV luminosities of the knots. Although Kraft et al. (2002) derived their LF for 0.1- 
10 keV luminosities of the knots with a different choice of spectral parameters (Fx = 2.3 and Nu = 
1.7x10^^ cm~^), a rough conversion to their 0.1-10 keV fiuxes can be obtained by multiplying our 
0.5-5 keV luminosities by a factor of three. 

The limiting sensitivity for this analysis is somewhat uncertain due to source confusion and 
possible variations of spatial extent among the jet-knots. Furthermore, the sensitivity of the ob- 
servation varies across the field of view due to the telescope vignetting and variations in the PSF. 
A detailed computation of the limiting sensitivity requires a Monte Carlo simulation taking into 
account the unknown morphology of the intrinsic diffuse emission and is beyond the scope of this 
paper. Here we follow the simplified method adopted by Kraft et al. (2001): we use the measured 
background and the telescope PSF at the edge of the field of view and we require a 4 cr measure- 
ment of the count rate. Our conservative estimate (using the background map and the telescope 
vignetting; Kraft et al. 2001) provides a sensitivity limit of Lgcus ~ 4x10'^^ erg s~^. Let us briefly 
discuss this value. During these observations the encircled energy of a point source is contained 
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within a diameter of ~3" all along the jet. Assuming that the intrinsic spatial extent of the knots is 
less than 4", we conservatively expect that all the compact sources can be detected within a circle 
of ~ 6" diameter, even if they are in the outermost part of the jet (i.e., > 4' off axis with respect 
to the mirror focus). The typical background within this region is ~15 cts from the observational 
data. To detect a knot at the required statistical level, we therefore need ~25 net source counts 
excluding the background photons. This corresponds to an X-ray luminosity of i^sens = 2x10^^ erg 
s^^. If a similar knot is embedded in extended emission (as observed in 50" < r < 80"), ~50 net 
photons are required over the ~120 background photons. It is therefore reasonable to assume -Lsens 
= 4x10^^ erg s~^ in the following discussion. 

In Figure 8, the hatched region shows where our sampling is incomplete and somewhat biased. 
The overall LF is well modeled by a broken power-law, with the power-law slope flatter below the 
break. The dashed line on Figure 8 shows the fit considered hereafter, where the power-law slope 
[k of N{> Lx) oc Lx**) flattens significantly from k = 1.4 to 0.5 below Lbrk — 9x10^^ erg s~^. 
The maximum likelihood method of Crawford, Jauncey & Murdoch (1970), using separate data 
sets below and above the break, supports this fit with an uncertainty (defined here as a standard 
deviation) of o" ~ 0.4 for each k. However, due to the limited number of the detected jet-knots (24 
for Lx > 4x10'^^ erg s^^), the Kolmogorov-Smirnov test (Press et al. 1992) indicates that a single 
power-law model with k~ 0.7 is also statistically acceptable. Even so, a single steep power-law with 
K > 1.0 can be rejected at the >90% confidence level. Assuming further that the distribution of 
knot luminosities extends down to some minimum value L^am, the integrated total X-ray luminosity 
can be evaluated quantitatively for both the resolved and the unresolved jet-knots (see Kraft et al. 
2002 for the detailed formalism). We note that the total luminosity of the Ccn A jet (i.e., the sum 
of the jet-knots and the unresolved diffuse emission) is 1.6x 10^^ erg s^^. Meanwhile, by integrating 
the LF down to any value of L^i^ we get < 8.1x10^^ erg s~^, at most, if we assume the "best-fit" 
broken power-law LF described above. Note that this is only 25% larger than the sum of the already 
resolved jet-knots. The situation is almost unchanged if we take a single power-law function with 
K ~ 0.7 for the LF. We therefore conclude that about half of the total jet emission, Lx ~ 8x 10^^ erg 
s~^, is truly diffuse in nature. 

We finally consider an apparent discrepancy between our results and previous reports. Kraft 
et al. (2002) suggested that the diffuse X-ray emission of the Cen A jet may be entirely explained 
by unresolved knots only if L^nin = 3x lO^^crgs^^ and there is no break in the LF below the 
sensitivity limit of Lgens ~3xl0'^^erg s~^ (corresponding to Lx ~lxl0^^erg s~^ in our LF). There 
are several possibilities to account for this discrepancy. First, the sensitivity limit adopted in Kraft 
et al. (2002) is calculated for a single ~35 ksec observation. As discussed in detailed in Kraft et al. 
(2001), the sensitivity is limited by the source luminosity, not by the background, over the most of 
the field of view. Therefore, our analysis based on the co-added ACIS-S image (95 ksec exposure in 
total) provides a factor of ~ \/3 deeper sensitivity for the source detection. Moreover, the PSF in 
the inner jet is a factor of ~2 smaller than previous observations (see §2), which further improves 
the signal-to-noise ratio. Second, by applying the same sensitivity limit as Kraft et al. (2002), it 
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would be difficult to see a flattening in the low-luminosity part of the LF. Thus, a steep power-law 
LF is consistent with the work of Kraft et al. (2002; who obtained k=1.1), although is implausible 
when the low-luminosity range (Lx < 10^'^ erg s~^) is taken into account. Third, it seems that 
quite different conclusions can be reached for different choices of k. In fact, integration of the LF 
down to arbitrarily small Lmin approaches infinity if ac > 1. Future deep observations of the Cen A 
jet are therefore strongly encouraged, although it seems to us reasonable that a significant fraction 
of the X-ray jet emission in this source is actually diffuse. 



4. Discussion 
4.1. The Jet 

The observed X-ray spectral index of the diffuse emission within the whole jet — including 
both the central flow and the limb-brightened edges — is steep, ax — 1.0-1.2. The implied energy 
distribution of the X-ray synchrotron emitting electrons therefore has to be Ueij) oc 7'"'' with 
s = 2 ax + 1 = 3.0-3.4, where 7 is the electron Lorentz factor. Such a spectrum is very similar to 
the one observed in the compact knots, which are believed to be the sites of shock waves, suggesting 
that the same, or a similar, particle acceleration process is acting in both regions. In the case of 
the shock acceleration, however, one expects the energy spectral index of the freshly accelerated 
electrons to be Sin = 2.0 in the strong non-relativistic regime, or Sin>2.0 when relativistic plasma 
velocities are considered (see, e.g.. Kirk & Duffy 1999, Ostrowski 2002). Thus, if first order Fermi 
shock acceleration is indeed responsible for formation of the X-ray (synchrotron) emitting electrons 
within the Cen A jet, the observed particle energy spectrum — both within and outside the resolved 
knots — has to be already modified by the (dominant) synchrotron losses, i.e., the spectral index 
of freshly accelerated high-energy electrons has to be increased by the radiative cooling |7|syn = 
(Tt 7^ -B^/Gtt meC to s = Sin + 1 (the 'strong cooling regime' with homogeneous magnetic field B; 
see Kardashev 1962). In fact, the appropriate cooling timescale is very short, tgyn = 7/|7|syn ~ 
20 sr^Q^^^ yrs (for a jet magnetic field i?_4 = i?/10~^G ~ 1, an observed synchrotron photon 

energy eio = e/lOkcV ~ 1, and a jet Doppler factor 6 ~ 1), consistent with this interpretation. We 
note in this context that the detailed analysis of the radio emission from a number of FR I sources 
of Young et al. (2005) revealed a 'universal' injection spectral index for the radiating electrons of 
Sin = 2.1, in accordance with the above discussion. We also note that synchrotron cooling within 
a non-uniform magnetic field can result in a spectral break larger than the standard one discussed 
above, s > Sin + 1 (e.g., Bicknell &; Begelman 1996). 

On the other hand, the very short synchrotron lifetime and the extended character of the 
diffuse X-ray jet component discussed here indicate that the appropriate acceleration process for 
the X-ray emitting electrons cannot be connected exclusively with the localized compact knots of 
the kpc-scalc jet, but instead has to operate effectively throughout the whole jet volume. Indeed, 
the maximum propagation length of the jet electrons emitting lOkeV synchrotron photons within 
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the 100 //G jet (equipartition) magnetic field can be estimated as Zrad = cFtrad ~ 10 pc for F ~ 2, 
which clearly indicates the need for the contimioTis acceleration of the radiating electrons in between 
the compact knots. Let us mention at this point, that the analysis of the expected 7-ray emission 
from kpc-scale FR I jets presented by Stawarz et al. (2005) suggests that sub-equipartition magnetic 
fields (which could increase Zj-ad) are not likely in these objects. The problem could disappear if 
many small scale strong shocks are present (below the fiux resolution limit for our point source 
detection), possibly generated at the supersonic and highly turbulent jet boundary layer. However, 
as argued in the previous section, based on the constructed luminosity function for the resolved 
jet-knots, the observed extended X-ray emission is most likely truly diffuse in nature, and not just 
the pile-up of the unresolved compact knots. 

A physical process which can naturally account for the distributed continuous acceleration of 
the radiating electrons within the whole jet volume consists of interaction of the particles with 
turbulent MHD waves. Indeed, by applying quasi-linear theory for particle-wave interactions (see, 
e.g., Schlickeiser 1989), one can find that with a reasonable set of physical parameters for the 
Cen A jet, the available maximum energies of the electrons accelerated in this way may be high 
enough to alow for efficient production of the synchrotron X-rays. On the other hand, the resulting 
electron energy spectrum is expected to be very flat: with a one-dimensional power spectrum of 
Alfvenic turbulnce W(k) (xk~'^ and a 'typical' range of turbulence spectral index q £ (1,2) — 
usually obtained in phenomenological models for the astrophysical MHD turbulence (see, e.g., the 
recent review by Muller &; Biskamp 2003, and references therein) — the steady-state solution for the 
momentum diffusion equation is a power-law with = q — 1 between and 1 (see, e.g., Borovsky 
& Eilek 1986 or Park & Petrosian 1995). However, the observed X-ray spectral index of the diffuse 
component of the Cen A jet is steep, requiring an electron energy slope of s ~ 3. In particular, 
taking the expected injection spectrum with Sin either or 1, the strong (synchrotron) cooling 
regime condition gives an observed electron energy index of s = 2 in both cases (see Aharonian 
2002), or an expected X-ray power- law slope ax = 0.5, i.e., much flatter than observed. It is 
important to note, however, that the quasi-linear theory for particle-wave interactions is of limited 
applicability for physical conditions within extragalactic jets, which may be highly turbulent, even 
with a possibly dynamically important pressure due to ultra-relativistic particles at the jet boundary 
(see Ostrowski 2000). In particular, the expected large ratio of the escape to the radiative losses 
time scales within the jet (see Stawarz Sz Ostrowski 2002) means that the stationary solution may 
not be realizable due to accelerated particle back-reaction, or that the turbulent diffusion could 
dominate the escape process. In such cases the quasi-linear approach may no longer be valid. In 
addition, radiative effects within a non-uniform magnetic field can modify the high-energy part of 
the electron energy distribution in a non-trivial way. A set of such problems is under investigation. 
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4.2. Hard Outer Sheath? 

Our analysis indicates the possible presence of hard X-ray features within the outer (fragmen- 
tary) sheath of the Cen A kpc-scale jet, extending beyond the X-ray limb-brightened jet edges. 
The X-ray spectral properties of the outer jet boundaries are most probably very different from 
the ones characterizing the central jet flow. In particular, the observed X-ray spectral index can 
be as flat as ax < 0.5. In addition, the physical connection of these regions to the main flow is 
not clear, although they do roughly coincide with the edges of the radio structure. One may note 
in this context that detailed studies of the extragalactic radio jets using the 'spectral tomogra- 
phy' technique has also revealed the presence of some features along the outer jet boundaries. In 
particular, Katz-Stone & Rudnick (1997) showed that the large-scale jet in the FR I radio galaxy 
3C449, for which the radio spectral index or ~ 0.5 is roughly constant along the whole jet, is 
surrounded by a steep-spectrum radio 'sheath'. Also, large-scale jets in the FR II radio galaxy 
Cygnus A are placed in the cdgc-brightcncd 'channels', which are consistent with the enhancement 
of radio-emitting electrons at the outer jet boundaries (Katz-Stone & Rudnick 1994). 

Unfortunately, the present X-ray data are not sufficient to convincingly demonstrate the pres- 
ence of a hard sheath around the Cen A jet. Let us however mention briefly a possible scenario which 
can account for formation of such hard features, namely synchrotron emission of the high-energy 
electrons which lose energy predominantly due to inverse-Compton radiation in the Klein-Nishina 
regime (hereafter 'IC/KN' process; see Aharonian &: Ambartsumyan 1985; Dermer & Atoyan 2002; 
Kusunose & Takahara 2005; Khangulian &: Aharonian 2005; Moderski et al. 2005). An interesting 
aspect of this scenario is that the observed electron spectrum can be flatter than the injected one, 
with the asymptotic power-law slope s « Sjn — 1- Let us therefore assume that the acceleration pro- 
cess similar to the one operating within the jet acts also within the outer boundaries of the outflow, 
where the magnetic field intensity is possibly lower than within the jet. Thereby the high-energy 
electrons cool partly due to the synchrotron emission, producing X-ray photons, but mostly due to 
the IC/KN radiation at 7-ray frequencies. The dominant IC/KN process thus flattens the injected 
electron spectrum, leading to the observed flat X-ray (synchrotron) continuum. 

7 —1/2 

The Lorentz factor of electrons emitting 0.1 keV synchrotron photons is 7x ~ 10 . Hence, 

the required energy of the seed photons to be upscattered by these electrons in the KN regime is 
£0 > meC^/47x ~ 10~^ b]/^ eV. At the same time, the energy density of the seed photons has to be 
C/q > t^B ~ 4 X 10~^° Bi4ergcm~^ in order to ensure the dominance of the IC/KN losses. Taking 
the lowest expected value S_4 ~ 0.01, this already excludes the cosmic microwave background 
radiation as the source of seed photons. However, Cen A is a giant elliptical galaxy known for 
its strong infrared emission: S'ir = 23 Jy at Air = 12 fim (Knapp et al. 1990). Assuming that all 
of this radiation is produced within r = Ikpc from the galactic center (see Israel 1998), one can 
obtain a rough estimate for the energy density of the galactic ~ 0.1 eV photons at the position of 
the kpc-scale Cen A jet, Uj^ = Smdf^/r^ Air ~ 2 x 10~^^ ergcm"^, close to the average value for 
giant ellipticals, ~ 10"^-^ ergcm~^, found by Stawarz et al. (2003). Note that this is about two 
orders of magnitude lower than the expected energy density of the Agtar ~ 1 /wm starlight energy 
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density (Stawarz et al. 2003). Also, the energy density of the synchrotron infrared emission of 
the nuclear jet in Cen A can be important in this respect. In particular, taking the observed 
synchrotron luminosity of the Cen A nuclear jet Lbi ~ 3 x 10^"*^ erg s^^ and its critical observed 
synchrotron frequency « 10^^ Hz (Chiaberge et al. 2000), one gets the appropriate energy density 
as measured by a stationary observer located on the jet axis at r = 1 kpc from the active center 
Uu = Lu rbi/47r c ~ lO"^*^ erg cm"^, and the appropriate critical frequency Fbi £bi ~ 0.03 eV, for 
the expected bulk Lorentz factor of the nuclear jet Fbi ~ 10 (see Stawarz et al. 2003). Relative 
enhancement of synchrotron emission produced within a structured relativistic jet may constitute 
an alternative/additional source of the seed infrared photons (see in this context Ghisellini et al. 
2005). Thus, we conclude that if only the magnetic field in the outer sheath of the kpc-scale Cen A 
jet is S 10 /xG can the suggested hard X-ray emission result from synchrotron radiation of the 
high energy electrons characterized by the standard initial spectrum Sin — 2, which cool mainly by 
Comptonization of the nuclear jet emission (and also the starlight photon field) in the KN regime.^ 



5. Conclusions 

Wc have presented an analysis of the diffuse X-ray emission from the kpc-scale jet in Cen- 
taurus A. We found that the extended emission is most likely diffuse in nature rather than being 
the summation of unresolved compact sources and that it exhibits a relatively flat-topped intensity 
profile with a steep, uniformly distributed spectral index of ax ^ Ij similar to those observed in 
the compact knots. We interpret the observed X-ray morphology in terms of a stratified jet model, 
consisting of a relativistic outflow with an expanding turbulent boundary shear layer. We argue 
that the observed X-ray spectrum of this diffuse component, consistent with the synchrotron emis- 
sion produced in the strong cooling regime, requires continuous acceleration of ultra-relativistic 
electrons within the whole jet volume. We also note that the jet may be (partially) surrounded 
by an outer sheath with a very hard X-ray spectrum. Both the presence and origin of this feature 
are unclear: we suggest that it forms in a turbulent region just outside the jet, separating the 
main outflow from the external medium, where the conditions for the radiative cooling processes 
are different from the ones characterizing the jet. Although the present data do not allow us to 
unequivocally demonstrate the presence of the hard outer sheath, or to determine the appropriate 
particle acceleration mechanism acting within (and possibly outside) the outflow, our analysis indi- 
cates that (1) the jet structure is more complex than that visible on the intensity maps alone, and 
(2) the acceleration and spectral evolution of the radiating particles may be governed by many dif- 
ferent processes (including those such as magnetic field reconnection which we have not discussed), 
linked to the complex magneto-hydrodynamical configuration of the outfiow. 



^We note that the very high energy 7-ray flux of the outer sheath impUed by this model, ~ 10 erg cm ^ s 
is weU below the sensitivity limits of current Cherenkov Telescopes. 
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Table 1: Newly dotcctcd X-ray knots and onhancomonts in the Con A jot 





Distance 


Angle 




R.A. 


Decl. 


Net Counts 


Knot ID 


(arcscc) 


(deg) 


(J2000) 


(J2000) 


(counts) 


1 


9.96 


36.65 


13 


25 


28.32 


-43 01 02.8 


159±18 


2 


15.31 


35.23 


13 


25 


28.73 


-43 00 59.9 


2682±55" 


3 


35.56 


36.96 


13 


25 


30.16 


-43 00 47.6 


244±17^ 


4 


40.96 


39.41 


13 


25 


30.45 


-43 00 43.1 


100±12 


5 


55.65 


35.06 


13 


25 


31.69 


-43 00 37.2 


204±19 


6 


65.96 


31.96 


13 


25 


32.63 


-43 00 34.3 


123±13 


7 


80.78 


40.48 


13 


25 


33.12 


-43 00 17.0 


42± 8 


8 


89.48 


40.05 


13 


25 


33.75 


-43 00 12.0 


51±10 


9 


86.64 


37.35 


13 


25 


33.79 


-43 00 16.9 


53±10 


10 


82.92 


32.40 


13 


25 


33.89 


-43 00 24.9 


44± 9 


11 


107.94 


31.10 


13 


25 


35.90 


-43 00 13.8 


46± 9 


12 


163.57 


38.42 


13 


25 


39.10 


-42 59 28.6 


16± 5 


13 


209.89 


37.59 


13 


25 


42.52 


-42 59 02.6 


38± 8 



"This extremely bright knot is only 1.5" from the neighboring bright knot (RA=13:25:28.57, Dec=-43:01:00.1) and 
therefore was not previously able to be resolved from AXl or knot2 in Table 2 of Kraft et al. 2003. 
''This jet-knot may be identified with the X-ray star J132530. 1-430048 reported in Colbert et al. (2004). 

Note. — Thirteen new X-ray jet-knots and enhancements have been detected, in addition to those in Table 2 of 

Kraft et al. (2002), in the co-added, 0.4— S.OkoV ACIS-S image. However, three knots in Kraft ot al. (2002), namely 
NXl (Knot-1), EX3 (Knot- 19), and FX4 (Knot-23), were not detected at their previously reported level for various 
reasons. NXl is only <3" apart from the nucleus and can hardly be detected in our 0.4—8.0 keV bandpass due to 
strong contamination of the bright nucleus omission. Kraft ct al. (2002) claimed a detection of NXl using only the 
0.4—1.5 keV bandpass image, since the heavily absorbed nucleus is expected to be suppressed. This is also the case 
in our analysis, but we adopted a more conservative approach. We find no evidence for any enhancement in the jet 
omission around EX3, suggesting it may have been a statistical fluctuation. Finally, wo suggest that FX4 is identical 
to FX3 reported in Kraft et al. (2002), but was regarded as a separate jet-knot due to poor photon statistics. 
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Fig. 1. — top: The 0.4-8.0 keV X-ray image of the Cen A jet. Archival data from Obs ID 2978 
and 3965 are combined. The image is smoothed with a a = 0.5 arcsec Gaussian, bottom: The 
0.4— 8.0keV X-ray image of the Cen A jet, after subtracting the point sources and fiUing in the 

blanks with values interpolated from surrounding background regions. More details are given in 
the text. Definitions of the angle and radial distance r are given in the panel. 
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Fig. 2. — Longitudinal jet profile of the Cen A in the soft X-ray energy band (0.4-1.5 keV; upper) 
and the hard X-ray energy band (1.5-8.0 keV; lower), where the photons counts are integrated over 
the position angles 25° < G < 45°. The dotted histogram shows the profile derived from the raw 
X-ray image (Figure 1 top), whereas the solid histogram with data points shows the one obtained 
for the diffuse emission only (Figure 1 bottom). The dashed line shows the background profile 
compiled from 15° < O < 25° and 45° < 6 < 55° regions. 
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Fig. 3. — Upper: The observed longitudinal X-ray intensity profile of the diffuse jet component in 
the 0.4— S.OkeV photon energy band. Photon counts are integrated over the position angle ranges 
25° < e < 45° (histogram), 25° < 9 < 33° (sohd hne), 33° < G < 37° (dashes), and 37° < 9 < 45° 
(hatched). Lower: Variation of the hardness ratio along the main jet axis. Hardness is defined 
as the ratio between the X-ray counts measured in the soft X-ray energy band (0.4— l.SkeV) and 
the hard X-ray energy band (1.5— 8.0 keV). Note that the hardness ratio changes significantly at 
the distance rc ^50", but stays almost constant at larger distances. The dotted line shows the 
hardness changes for the background compiled from 15° < 9 < 25° and 45° < 9 < 55° regions. 
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Fig. 4. — Variation of the X-ray flux density measured at IkeV (top), the X-ray energy spectral 
index ax {middle), and the absorbing column density A^h {bottom) of the diffuse emission along 
the jet. Note that Nu changes significantly at Tc ~50", but that the energy spectral index is 
almost constant (ax — 1) along the jet. Hatched regions show the best-fit parameters and their 
la uncertainties for the combined spectral fittings within 10" < r < 50" and 50" < r <250", 
respectively. 
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Fig. 5. — Transverse jet profile of the Cen A in the soft X-ray energy band (0.4-1.5 keV) and the hard 
X-ray energy band (1.5-8.0 keV). The profiles for the inner jet region (upper, (a) 50" < r < 170") 
and the outer jet regions {lower; (b) 170" < r < 250") are given in separate panels. The dashed 
histogram shows the profile derived from the raw X-ray image (Figure 1 upper), whereas the solid 
histogram with the data points shows that obtained for the diffuse emission only (Figure 1 bottom) . 
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Fig. 6. — Upper: The observed transverse X-ray intensity profile of the diffuse jet component in the 
0.4—8.0 keV photon energy band. Lower: Variation of the hardness ratio in the direction perpen- 
dicular to the jet main axis. Hardness is defined as the ratio between the X-ray counts measured 
in the soft X-ray energy band (0.4— 1.5keV) and the hard X-ray energy band (1.5— 8.0 kcV). The 
profiles are given in separate panels for various distances from the nucleus, (a) 50" < r <80", 
(b) 80" < r <120", (c) 120" < r < 170", and (d) 170" < r <250". Note the possible hardness 
variations near the edges of the jet, 6 ~ 25° and 9 ~ 45°. 
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Fig. 7. — Variation of the X-ray intensity measured at IkeV (top), and the energy spectral index 
{bottom) across the jet, in A© = 2° steps. The contribution of point sources is simply removed, 
and dmfilth is not used in order to avoid any possible artifacts in the spectral analysis. Reduction 
of integration regions due to the point source subtraction is corrected in estimating the 1 keV flux. 
The fitting model consists of a power-law function and the assumed Galactic absorption column 
density Nn = 0.96 x lO^^ cm'^. 
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Fig. 8. — Luminosity function of the knots in the jet, constructed excluding the knots of the 
innermost part of the outflow (0" < r <50"). The X-ray luminosity is the absorption corrected 
luminosity in the 0.5— 5.0keV band, assuming a power-law spectrum with energy spectral index 
ax = 1-0 and A^h = 0.96x10^^ cm~^. The best fit broken power-law model, excluding the hatched 
region where our sample is incomplete and biased (Lx <4xl0^^erg s~^), is also shown. 



